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Abstract 



The paper "Big Bang Nucleosynthesis and Active-Sterile Neutrino Mixing: 
Evidence for Maximal y„ *-* v T Mixing in Super Kamiokande?" by X. Shi 



and G. M. Fuller ( astro-ph /9810075 ) discusses the cosmological implications 
of the Vu — ► v s solution to the atmospheric neutrino anomaly. It incorrectly 
concludes that a lower bound on the v T mass of 15 eV is needed in order for the 
v a —> v s oscillations to be consistent with a Big Bang Nucleosynthesis upper 
bound on the effective number of neutrino flavours of 3.3. Since such a large 
v T mass is disfavoured from large scale structure formation considerations, the 
strong, but incorrect, conclusion is made that cosmology favours the — > v T 
solution to the atmospheric neutrino problem. We explain the nature of the 
error. We conclude that cosmology is, at present, consistent with both the 
Vn, — ► v s and Vn v T possibilities. 
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The large up-down asymmetry observed by SuperKamiokande for atmospheric neutrino 
induced /i-like events provides compelling evidence for the disappearence of muon-neutrinos 
The most natural explanation of this phenomenon is neutrino mass driven oscillations 
between and another flavour u x , with oscillation parameters in the approximate range 

1(T 3 < AmJJeV 2 < 10" 2 , 0.8 < sin 2 29^ < 1, (1) 

where Am 2 ^, is the squared mass difference, and 9^ x is the vacuum mixing angle. The 
SuperKamiokande results cannot at present distinguish between the competing v x = v s and 
v x = v T possibilities, where v s is a sterile neutrino. 

The v s solution appears to provide a cosmological challenge: with the parameters as 
per Eq.(|l]) naive calculations |§ suggest that the v s is thermally equilibrated prior to the 
Big Bang Nucleosynthesis (BBN) epoch, leading to an effective number of neutrino flavours 
iVeff = 4. It is at present unclear if N e Q = 4 is inconsistent, but the low primordial deuterium 
abundance result reported in Ref. favours N e g <C 4. The naive calculations referred to 
above ignore the creation of neutrino-antineutrino asymmetries by active-sterile neutrino 
oscillations 0. Large neutrino-antineutrino asymmetries suppress z/ M — > v s oscillations, and 
generically invalidate conclusions drawn from the naive calculations. 

Neutrino asymmetry evolution as driven by active-sterile oscillations goes through two 
distinct phases at temperatures immediately below the critical temperature T c at which 
neutrino asymmetry growth begins. The first phase is explosive growth that initially is 
exponential in character. After a time, the non-linear evolution equation for the neutrino 
asymmetry turns the exponential growth into power law T -4 growth. 

We have shown in Ref. that small angle v T — > v s oscillations with m Vr > m Us will, 
for a range of parameters, create a large v T asymmetry which will then suppress — > u s 
oscillations. The situation is a little complicated because — > v s oscillations act to create a 

asymmetry which can compensate for the effect of the large v T asymmetry in the effective 
potential for — + v s oscillations. Indeed the Wolfenstein part of the effective potentials 
|7j for v T — > v s oscillations and — > v s oscillations are respectively proportional to the 
quantities 

~2L VT + L Vll , 

LM~2L Vii +L Vt , (2) 

where Lf = (nf — nj)/n 7 is the asymmetry for species / whose number density is denoted 
by nf. The generation of L Vt asymmetry by v T — > v s oscillations will also generate a large 
asymmetry, which then suppresses — > v s oscillations provided that these oscillations 
do not generate a significant L v asymmetry in such a way as to drive to zero. Whether 
or not the v s oscillations can compensate for the effects of the large v T asymmetry will 
depend on the rates at which L Ut and L u are generated. The rates at which L Ut and L u 
are generated depend on the oscillation parameters, Am 2 and sin 2 29, for each oscillation 
mode. If we fix Am 2 s and sin 2 29^ s to fit the atmospheric neutrino anomaly then this leaves 
two free parameters, Am 2 s and sin 2 29 TS for the v T — * u s oscillations. The issue therefore 
is to calculate the region of (sin 2 29 TS , Am 2 J parameter space where the large v T — > u s 
oscillation generated v T asymmetry is not compensated by a — * v s oscillation generated 
asymmetry. For reasons explained in a moment, this issue is best handled numerically. 
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The numerical integration of the coupled evolution equations for L Vt and L Vfl for this system 
was first done in Ref. |]], and updated in Ref. ||. See Figure 7 of Ref. || and Figure 2 
of Ref. || for the pertinent results. These results show that m Ur > few eV is required, 
with the precise value being a function of sin 2 2# TS . Thus, the interesting conclusion is that 
if the atmospheric neutrino anomaly is due to —>■ v s oscillations then the v T must have 
a cosmologically interesting mass. Furthermore, the mass of the v T should be in the range 
which can be probed by the Chorus and Nomad experiments. 

Also contained in Ref. |J (see section VI) was an approximate analytic computation of the 
allowed region of parameter space. The main difficulty in solving this problem analytically 
is that the evolution equations for L Vt and L„ are complicated coupled non-linear integro- 
differential equations. For the approximate analytic computation of Ref. || we made the 
simplifing assumption that the v T asymmetry is in the T -4 phase when — > v s oscillations 
are most strongly acting to produce a compensating L v . This simplifing assumption makes 
things much easier since it means that 

dL Ur -AL Ut 

With this simplifying assumption (plus some other assumptions) we obtained the allowed re- 
gion of parameter space where the v T — > v s oscillation generated L Vt asymmetry successfully 
suppressed maximal — > v a oscillations. The result was the following simple expression 
[see Eq.(144) of Ref . gj] Q, 
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|Am 2 s /eV 2 | ^ 6 x 10 5 (|Am 2 s /eV 2 |) ' (4) 



This result is independent of sin 2 29 TS just because Eq.(|3]) is independent of sin 2 29 TS . How- 
ever, our numerical work (also in section VI of Ref. ||) showed that the simplifying assump- 
tion that leads to Eq.(^) is not valid. The reason is that the v T asymmetry is still in the 
explosive growth phase when — > v s oscillations are most strongly acting to produce a 
compensating L u asymmetry. This means that the rate at which L Vt is generated is actu- 
ally much greater than assumed in the derivation of Eq.(f|). This makes the allowed region 
significantly larger and it also depends on sin 2 29 TS (since the rate of L Ut generation depends 
on sin 2 29 TS in the explosive growth phase). Thus, it turns out that the analytic approach 
is not very accurate because it underestimates the allowed region by up to several orders of 
magnitude in Am 2 s . However, it does provide useful qualitative insight, which is why we 
included it in Ref. 0. 

This analytic approach has recently been re-examined by Shi and Fuller in Ref. || . Using 
a similar set of approximations, including Eq.(|^), they find the slightly less stringent bound 
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lAmyeV^SxlO^IAm^/eV 2 !)' 7 . (5) 



Note that the meaning of this bound, assuming for the moment that it is exactly valid, is as 
follows: If it is obeyed, then negligible numbers of sterile neutrinos are produced by — > v s 
oscillations. If it is not obeyed, then the sterile neutrinos will eventually be brought into thermal 
equilibrium by — > v s oscillations. 
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Taking Am^ s = 10~ 3 eV 2 , Eq.(||) then implies a lower bound on m Vr of about 15 eV. 
However, as we have explained above, and as was also pointed out in Ref. |J (see pages 5167 
and 5170), this analytic approach significantly underestimates the allowed region because 
the approximations upon which it is based are simply not valid. The main problem, as 
discussed above, is the assumption that the v T asymmetry is in the T -4 phase when — > u s 
oscillations are most strongly acting to produce a compensating L v asymmetry. 

This point is easy to demonstrate. Consider, by way of concrete example, v T — > u s 
oscillations with the parameters Am^ s = —50 eV 2 and sin 2 29 TS = 10~ 8 . In the absence 
of other oscillation modes involving either v T or u s , the growth of the v T asymmetry is 
illustrated in Fig.l of Ref. || which we reproduce here for convenience. In this example the 
critical temperature is about 37 MeV, and the value of the tau-like asymmetry L Vr is about 
10~ 5 when the explosive growth phase gives way to T -4 growth. 

Consider now maximally-mixed — > v s oscillations. The MSW resonance condition for 
this mode is given by 

V(p) = Vwolfenstein + ^finite-T = (6) 

where 

Vwolfenstein - ±V%G F n 7 (2L Ufl + L„ T ) (7) 

and 



s/2G F n 7 A(—) (8) 
\m w / (v) 



finite-T 



,m w J (p) 



are the Wolfenstein (finite-density) and finite-temperature contributions to the effective mat- 
ter potential 0, respectively [the + (— ) sign refers to neutrino (anti-neutrino) oscillations]. 
The quantity G F is the Fermi constant, A ~ 15.3, my/ is the IV-boson mass, p is the neutrino 
momentum, and (p) ~ 3.15T is its thermal average. Taking for definiteness that L Ut > 0, 
which means that the MSW resonance for — > v s oscillations occurs for the neutrinos 
rather than the anti- neutrinos, Eq.(^) can be solved for the resonance momentum p rcs to 
yield 

(p) ~ AT 2 K ' 

provided that L V)1 L Vt . Putting in Lp T ~ 10 ^ and T ~ 37 MeV we find that 

Pres Pre; 



10. (10) 



This means that as L Vt develops from close to to about 10~ 5 during the explosive growth 
phase, the MSW resonance momentum for — > v a oscillation sweeps through most of the 
Fermi-Dirac distribution. This is what we mean when we say that the tau-like asymmetry 
is still in the explosive growth phase when — > v s oscillations are most strongly acting to 
produce a compensating L v asymmetry. 
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Of course, the above was derived assuming that — > u s oscillation are not strong enough 
to create a significant compensating L v , so that L v <C L Vt always holds. This is of course 
the main point. Are — > v s oscillations effective enough so that they create significant L v 
so that p ms ~ (p) contrary to Eq.(|TU|)? In order to work out exactly what happens, a detailed 
numerical calculation must be performed, which is one of the points we were emphasising 
above. However, the above argument clearly shows that L Vt should not be taken to be in the 
T~ 4 phase. Note that if sin 2 29 TS > 10~ 8 then the explosive growth phase typically continues 
until even larger values of L Ur are reached (i.e. larger than 10~ 5 ). This obviously makes the 
T" 4 growth assumption even less valid. 

To summarise, — ► u s oscillations can solve the atmosheric neutrino problem without 
leading to N e g = 4 neutrinos in the early Universe. This result occurs because v T — > v s 
oscillations generate a large L Vr asymmetry (provided m„ r > m Va ). The large L Ut asymmetry 
will suppress the — > v s oscillations in the early Universe provided that the — > z/ s 
oscillations do not generate a large which compensates for the effect of the large L Vt in 
the matter term for — > v s oscillations. Calculating the 'allowed region' is best handled 
numerically and has already been done in Refs. P,j8|- The analytic approach of Shi and 
Fuller H, which is based on our own analytic work, significantly underestimates the allowed 
region because it uses simplifying assumptions which are not valid. We conclude that there 
is at present no cosmological objection to the —>■ u s solution to the atmospheric neutrino 
problem. 
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Figure caption 

Figure 1: The evolution of \L Ut \ for the example, Am^ s = —50 eV 2 , sin 2 29 TS = 10 
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